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Abstract —Recently, spatial stochastic models based on deter- 
minantal point processes (DPP) are studied as promising models 
for analysis of cellular wireless networks. Indeed, the DPPs can 
express the repulsive nature of the macro base station (BS) 
configuration observed in a real cellular network and have 
many desirable mathematical properties to analyze the network 
performance. However, almost all the prior works on the DPP 
based models assume the Rayleigh fading while the spatial models 
based on Poisson point processes have been developed to allow 
arbitrary distributions of fading/shadowing propagation effects. 
In order for the DPP based model to be more promising, it is 
essential to extend it to allow non-Rayleigh propagation effects. 
In the present paper, we propose the downlink cellular network 
model where the BSs are deployed according to the Ginibre 
point process, which is one of the main examples of the DPPs, 
over Nakagami-m fading. For the proposed model, we derive a 
numerically computable form of the coverage probability and 
reveal some properties of it numerically and theoretically. 

I. Introduction 

Spatial stochastic models for wireless communication net¬ 
works have attracted much attention (see, e.g., m-E). This 
is because the performance of a wireless network depends 
on the spatial configuration of wireless nodes and spatial 
point processes can capture the irregularity of the node con¬ 
figuration in a real network. In the analysis of such spatial 
models for wireless networks, the stationary Poisson point 
processes (PPPs) have widely used to model the configuration 
of wireless nodes. While such a PPP based model is tractable 
to analyze, it is an idealized one and ignores the spatial 
correlation among the wireless nodes. 

Recently, spatial models based on determinantal point 
processes (DPPs) are studied as promising models for analysis 
of cellular wireless networks (see, e.g., 0 -hsd- Indeed, the 
DPPs can express the repulsive nature among the macro base 
stations (BSs) observed in a real cellular network and have 
many desirable mathematical properties to analyze the network 
performance (see, e.g., IITTl - liT3l ). However, almost all the 
prior works on the DPP based models assume the Rayleigh 
fading while the PPP based models have been developed to 
allow arbitrary distributions of fading/shadowing propagation 
effects (see, e.g., m-m. There are only a few exceptions 
as follows: Although 0 mainly assume the Rayleigh fading, 
it states in Remark 2 that the fading from the interfering 


BSs can be generalized to follow an arbitrary distribution 
while retaining the Rayleigh fading from the serving BS. 
Also, 0 investigates the tail asymptotics of the distribution 
of interference using different fading distributions, though it 
does not take into account the fading from the serving BS. 

In order for the DPP based model to be more promising, 
it is essential to extend it to allow non-Rayleigh propagation 
effects from both the serving and interfering BSs. In the 
present paper, we propose the downlink cellular network model 
where the BSs are deployed according to the Ginibre point 
process (GPP) over Nakagami-m fading (see lfl8l ). The GPP 
is one of the main examples of the DPPs and is suited to model 
the macro BS configuration (see e.g., 0). The Nakagami- 
m fading includes the Rayleigh one as a special case of 
m = 1 and makes the fading effects to follow the mth 
Erlang distribution with unit mean. It is further important 
to study the Erlang distributed propagation effects since it 
could lead to an extension to the analysis of multi-input multi¬ 
output heterogeneous cellular networks as in fT9l — ll2ll . For the 
proposed model, we derive a numerically computable form 
of the coverage probability and reveal some properties of it 
numerically and theoretically. 

The paper is organized as follows: In the next section, 
the network model which we consider is defined, where the 
impact of the variability of fading effects on the signal-to- 
interference ratio (SIR) is discussed using some stochastic 
order. A brief description of the GPP is also provided. In 
Section [III] we first give the basic formula of the downlink 
coverage probability when the BSs are deployed according 
to a general stationary point process and then we derive a 
numerically computable form of it when the BSs are deployed 
according to the GPP. In Section [IV] the asymptotic properties 
of the downlink coverage probability in two extreme cases 
m = 1 and m —> oo are discussed. The results of some 
numerical experiments are shown in Section [V] 

II. Network model 
A. Macro base station network 

The downlink cellular network model which we consider 
mainly follows the ones in 0, ED and consists of homoge¬ 
neous macro BSs. Let $ denote a point process on M 2 and 


let Xi, i £ N, denote the points of <f>, where the order of 
Xi, X 2 , ..., is arbitrary. The point process $ expresses the 
configuration of the BSs and we refer to the BS located at Xi 
as BS i. We assume that $ is simple and locally finite a.s. 
and also stationary. Each mobile user is associated with the 
closest BS. Thus, due to the stationarity of the point process 
and the homogeneity of the BSs, we can focus on a typical 
user located at the origin o = (0,0). For the simplicity, we 
limit ourselves to the interference limited (noise-free) case 
throughout the paper. It is known that, in the interference 
limited case of the stationary single tier network model, 
the downlink coverage probability does not depend on the 
intensity of the point process, transmission power and path- 
loss coefficient (see, e.g., S3, ED), so that, we set these 
values as 7 r _1 , 1 and 1, respectively, for convenience. The 
path-loss function £ representing the attenuation of signals 
with distance is then given by £{r) = r~ 2 r > 0, for some 
f3 > 1. We assume the Nakagami-m fading, so that the random 
propagation effect from the BS i to the typical user, denoted 
by Hi, follows the ?nth Erlang distribution with unit mean, 
denoted by Hi ~ Gam(?n, 1/m), i £ N, where Hi, t £ N, 
are mutually independent and also independent of the point 
process <£. The shadowing is ignored in our model. 


In this setting, the downlink SIR of the typical user is given 

SIR. - « (1, 

where B a denotes the index of the BS associated with the 
typical user at the origin; that is, {B a = i} = {||JQ|| < 
|Xj||, j £ N}, and I 0 (i), i £ N, denotes the cumulative 
interference received by the typical user from all the BSs 
except BS i; that is, 

Io (*)= E H ^(\\ X s\\)- (2) 

When we emphasize the m of the Nakagami-m fading, we put 
the superscript (m) and write SIR/ m ), l/ m )(j) and so on. 


Before assigning a specific point process to the BS con¬ 
figuration <1>, we here discuss the variability of the SIR in 
terms of the value of m of the Nakagami-m fading using 
some stochastic order. For two random variables X and Y with 
finite expectations, we say that X is less than Y in the convex 
order [resp. increasing convex order] and write X < cx Y 
[resp. X <i cx Y\ if E f(X) < E f(Y) for any convex [resp. 
increasing and convex] function / such that the expectations 
exist (see, e.g., Il23i ). It is well known that Gam(m, 1/m) is 
decreasing in m in the convex order; that is, mi < m 2 implies 
f?( m 1 ) > cx for ~ Gam(m, 1/m). However, we 

have to care about the application of the convex or increasing 
convex order to the SIR in ([TJ since it may not have the 
finite expectation when the path-loss function is unbounded 
at the origin. We thus consider the conditional SIR given the 
common BS configuration <t>. The following lemma implies 
that more variable fading effects lead to larger and more 
variable conditional SIR given the common BS configuration. 

Lemma 1: For each k = 1, 2, let {///*") }ieN denote a 
sequence of mutually independent random variables represent¬ 
ing the fading effects. Let also SIR/ fc ), k = 1, 2, denote 
the SIR of the typical user, given by ©, corresponding to 


the fading effects {H- I> } ie m with the common BS configu¬ 
ration $ = {JGjigft- Then, H 1 -' 1 > cx h\ 2 \ i £ N, implies 
(SIR/ 1 ) I $) >icx (SIR/ 2) | $) given <1>. 

Proof: By Strassen’s theorem (see, e.g., (23] Sec. 1.5]), we 
can define {ir/ fe )} ieN , k = 1, 2, satisfying H l j ' 2 ' 1 = E(fT,- 
JT/ 2 )) a.s. for each i£N. Thus, by ([TJ and Jensen’s inequality, 

S|R (2) = y E(Hf) \H[ 2) )£(\\Xf\) 

£ j 6 N\ {i } \H™)t(\\x j \\) 

x 1 {ll x «ll<ll A 'jll, few} 

< E(SIR/ 1} I {Jrf)} ieN ,$) a.s. 

Let / denote any increasing and convex function satisfying 
E(/(SIR/ 1} ) I $) < 00 and E(/(SIR/ 2 )) | $) > -00 a.s. 
Then, Jensen’s inequality implies 

e(/(sir/ 2 )) 1 $) < e(/(e(sir/ 1 ) I I d>) 

< E(/(SIR/ 1) ) I $) a.s. 


B. The Ginibre point process 

We here give a brief description of the GPP (see, e.g., 
111 1 — B1 31 for details). The GPP is one of the DPPs on the 
complex plane C defined as follows. Let <I> denote a simple 
point process on C, and let p n : C" —> K.+, n £ N, denote 
its ?rth joint intensities (correlation functions) with respect to 
some Radon measure v on C; that is, for any symmetric and 
continuous function / on C n with compact support, 

e( E f(Xi,X 2 ,...,X n ) 

distinct 

I f{z\,Z 2 ,...,Z n ) 

C" 

x p n {zi,z 2 ,... ,z n ) v(dzi) v(dz 2 ) ■ ■ ■ v(dz n ). (3) 

The point process $ is said to be a DPP with kernel K\ C 2 —>• 
C with respect to the reference measure v if p n in (|3j is given 
by 

Pn{zi,z 2 ,... ,z„) = det (K(zi, Zj))?j =1 , 

where det denotes the determinant. Furthermore, a DPP $ is 
said to be the GPP when the kernel is given by K ( z,w ) = e zw , 
z, w £ C, with respect to the Gaussian measure is(dz) = 
7 r” 1 p(dz), where w denotes the complex conjugate of 

w £ C and p denotes the Lebesgue measure on C. It is known 
that the GPP is motion invariant (stationary and isotropic) and, 
by the definition, its intensity is equal to 7r -1 ; that is, E<f>(C) = 
7T _1 p(C) for C £ £>(C). One of the useful properties of the 
GPP is given by the following proposition. 

Proposition 1 (Kostlan J24K ): Let X tl i £ N, denote the 
points of the GPP. Then, the set {||X,|| 2 }i g N has the same 
distribution as where Y t , i £ N, are mutually 

independent and each Y t follows the vth Erlang distribution 
with unit-rate parameter, denoted by Y, ~ Gam(i, 1), i g N. 





III. Downlink coverage probability 


B. Ginibre base station deployments 


The downlink coverage probability is defined as the prob¬ 
ability with which the SIR of the typical user achieves a target 
threshold; that is, P(SIR 0 > 9) for 9 > 0. In this section, we 
first derive the basic formula concerning the downlink coverage 
probability when the BSs are deployed according to a general 
stationary point process, and then we provide its numerically 
computable form when the BSs are deployed according to the 
GPP. 

A. General stationary base station deployments 

Lemma 2: Consider the downlink cellular network model 
such that the BSs have mutually independent Nakagami-m 
fading channels and are deployed according to a stationary 
point process $ on R 2 . Then, the coverage probability for a 
typical user satisfies 


p(siRr > o) 

= E^ e 

• ^ n ! 


n—0 


d" 

dx n 


TT (i i a m 

n \ i+e mx^T) x ) 


jen\{B a } 


X = l- 

(4) 


Proof: Throughout the proof, we fix the value of to and 
drop off the superscripts (m). We have from 0 that 

P(SIR 0 > 9) = P(SIR 0 >9, B a = i) 

ie N 

91 0 {i) 


E P (E> 


ieN 


*(ll*ill) ! 


*11 < Ill'll. J6N). (5) 


For each i £ N, since Hi ~ Gam( to, 1/to), conditioning on 
$ = { X /}.?eN and {Hj} jen \ {i} yields 


'( 


PI Hi > 


91 0 (i) 

*(ll*ll) 


4 1 , {Hj}jen\{i} 


^ 1 rm9I 0 (i)\ n t m9I 0 (i ) 
^ n! W(M) ) “"l /(KID 


n—0 
m —1 


= E 


n—0 


(-l) n d n 
n\ dx 7 


■ exp 


(- 


*(ll*ll) 

m91 0 (i) 


*(ll*ll) 


‘) 


( 6 ) 


Here, let us take the conditional expectation given $. Since 
exp (—m91 0 (i) ir/£(||A,;||)) is continuous in x and each Hj ~ 
Gam(TO, 1 /to) has the finite moment of any order, we can 
exchange the conditional expectation and derivative. Thus, 
since Hj , j £ N, are mutually independent, applying © yields 


(s? exp (- 


_d 

^da; 

E 

d.T 


1 9I a (i) 




<(11*11) 

n E(exof 11 *11) 

«**»> 

jL TT 

d.T 11 


$ 


■ n (i+- 

jeN\{i} 


*(ll*ll) 


(7) 


of 


where the Laplace transform Ch(s) = (1 + s/m)~ 

Hj ~ Gam (to, 1/to) is applied in the second equality. Finally, 
applying (J6]» and 0 LO ©. we obtain 0. ■ 


We consider the GPP as the BS configuration ( l», where a 
point z = a; + iy£:Cis identified as (x, y) £ R 2 . 

Theorem 1: Consider the downlink cellular network model 
with path-loss function £(r) = r -2 4, r > 0, such that the 
BSs have mutually independent Nakagami-m fading channels 
and are deployed according to the GPP. Then, the coverage 
probability for a typical user is given by 

P(SIR^ m) > 9) 

poo m— 1 n r k 

= / e~ u #)(«) £ 9 n E E(-E M S £h(' 


n=0 k—0 L h=l 

n—k ^ rr—1 


(hi,h,2, 


e nj, D-vrf: 

■,h n -k)eVn-k r=1 r L 9=° 


( 3,(«) 


drt, 

( 8 ) 

where a h = a o- Pk = {(hi, h 2 , ■ ■ ■, h k ) £ (NU{0}) fe | 

Ej=i r h r = k}, and 


M^(u) = 

2=0 
CO l 

sfM = E TT (4"' 0) wr*' 1 rfilw, 

i=0 l - 

OO 


(9) 

GO) 

( 11 ) 


2=0 


with 


k-h-\-1 


«,<»)= e n / 

[ri,...,rk-h+l)&Qk,h 9—1 


T ^ 
' 9- 


/TO + q — 1\ 
\ to — 1 / 


[m,q) 


(u) 


( 12 ) 


Qk,h = {(n,r 2 ,...,r k - h+1 ) £ (N U {0}) fe_/i+1 | 

spk-h+l spk-h+1 .1 . 

Z^tq—l C 1 T q — K 5 2^/q—l r q ~ rL } * dnu 


/<"■»>(«)= 


e-‘v i Wv) ,f , 

-^dy, 9 = 0,1,2,.... 


(l + 9 (u/yY) 


(13) 

Remark 1: The numerical computation of (f8l> (IT3l> seems 
complicated. Note, however, that the infinite product and 
infinite sums in (Mi (fTTl i are not nested. The computation is 
thus scalable in the sense that these infinite product and infinite 
sums can be computed simultaneously in one iteration scheme. 
The case of to = 1 in Theorem Q] of course, coincides with 
Theorem 1 in 0. The cases of to = 2 and 3 are provided 
after the proof. 

Proof: As in the proof of Lemma [2j we drop off the 
superscripts (to) and write, for example, for j{ m ' q K By 
Proposition Q] we can set {||Xi||}j g N satisfying ||A^|| 2 ~ 
Yi ~ Gam(i, 1) for each i £ N, where Y), i £ N, are 
mutually independent. Thus, applying the density functions of 
Gam(i, 1), i £ N, and i(r) = r _2/3 , r > 0 to 0 yields 

P(SIR 0 > 9) 



























(17) 


m —1 


E^E^ n ■#+»(£)'.) 

ra=0 2=1 L jeN\{i} 


Yi\P \~* 
L 


X 1 


{E>Vi} 


K 


£c=l- 


00 -1 />oo 


_ (- 1 )" y, _l. 

n! i! 

n =0 2=0 


n 7 4 

3 ^n c >fe“) 


3=0 


dw, 


X =1 


(14) 


where 


Cj{x;u) = — 


e -v y j 


■ d y. 


j'Ju (l + 9(u/y)P x)' 

Note here that C :! (1; u) = Jj°\u) in (THT i. The general Leibniz 
rule (see, e.g., Il25l Sec. 5.2]) leads to 


1 f\j l Yh 

— IlQ(^w) = ^\{Ci(x;u)y 1 l[C j (x-,u) 


3=0 


daT 


£0 (i5) 


, -i • 
3=0 

i d n_fe °° 


k =0 


J=0 


Consider the first derivative on the right-hand side above. It is 
certainly ( Ci(x ; u )) when k = 0. For k = 1,2,..., Faa di 
Bruno’s formula (see, e.g., |j26l ) leads to 


-j^(<3(*;u)) 1 = Y J {~0 h h\{C i {x ] u)) 


-h -1 


h= 1 


X(*;«), ^ 2 ) (a:;u),...,C^ fc /t+ 1 ) (x;-u)), (16) 
where Bk,h denotes the Bell polynomial; 


Bk,h{p^l i ? • • • 5 /i+l) 

fc—Ji+1 1 

=‘ ! e n A®*. 


(ri,...,rk-h+i)&Qk,h 9=1 9 ^ 

and for 3 = 1 , 2 ,..., 


Cf 9) (x;w) = (-1) 9 0 9 
1 z* 00 


(m + 3 - 1 )! 

(to — 1)! 

e~y y 1 (u/y)^ 

m+q 


i- Ju (l + 9 ( u/y)P x ) 


dy. 


Note here that 
1 

9 




so that, noting that ^2 q=1 qr q = k and applying (fl 2 l> . we 
have 

B Kh (cf > (1;«), Cf } (1;, cf" fc+1) (1;«)) 

= (-l) fc 0 fc fc!E fc , M ( u ). 

Hence, applying this to ( IT 6 ] >. multiplying u l /i\ and summing 
over i = 0, 1, ... (see (fl4l >). we have 

°° n ,i ^/c 


Sirs vd c ^y' 


2=0 


= {-i) k e k k\Y J {~0 h h\s k , h {u). 


h=l 


On the other hand, we apply Faa di Bruno’s formula twice 
to the second derivative on the right-hand side of (IT5l) : first to 

I! C j ( x ;u) = exp (^ log Cj (x: u)) and then to log Cj(x;u). 

We then have 


i—k 


dx‘ 


IZfc n Cj{x\u) = (n — k)\ n Cj{x; r i 


3=0 


3=0 


i — k 


(hi,h 2 ,...,h n _i c )£V n -k r—1 


l h r 


E IIo hE^ 1osC7 ^ ; 


r! ^ dx r 
3=0 


(18) 


and 


— log Cj(x; u) = y](-l ) 9 1 (9-1)! (Cj(x-,u)y 


9=1 


x B r, q {cf\x\u), cf\x\u),...,C^ q+1 \x-,u)). (19) 
Here, taking x = 1 reduces to 

B r , q {cf ] (1; u), Cf\ 1; c\ s ~ r+1 \ 1; u)) 

= (-1 ) r r\e r V r , qJ (u), 

so that, applying this to ( fi~9] i and (fl~ 8 l > with Ylr=i r = n—k, 
we have 


jn-fe °° 

^ II C 3 (x; u) 


dx 


3=0 


= (-l)"- fc 9 n ~ k (n - k)\M(u) 


n—k ^ r r—1 


h r \ 

(hi,h 2 , — ,h„-k)eP n -k T '=l 9=0 


E lira El-OV^+.H 


( 20 ) 


Finally, we obtain © from (1 1 41 . (IT5] >. (IT7l) and (l 20 l >. ■ 

Example 1: We here give the examples of Theorem [T] for 
to = 2 and 3. When to = 2, it reduces to 

pOO 

P(SIR£ 2) > 9) = / e~ u M {2 \u) 

Jo 

x (45( u ) + 9 (45 ( u ) t i 2 i ( u ) - 45(«))) du - 


where 


aS 3 w = E: 
sgw = 2 E 


=s< ■“>(«)’ 

J ,~ ' O’M 
i!(jf' 0 , (o))' 


oo r (2,l)/ , 

«w= ! e 4 J- 

2—0 Jl m 

Also, when to = 3, it reduces to 

pOO 

P(SIRi 3) > 0) = / e~ u M {3) (u) 

Jo 

x sgj («)+ 9 (45 («) 45 («) - 45 («)) 
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where 



- sg (u) rg («) - sg («) + 2 sg ( U ) 


du, 
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^3(«) 

S'g (u) 

5g(«) 

TffM 

?fiV) 

T$(u) 
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oo 

3 E- 
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d./ 
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(m) 
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v 2 ’ 


6 E- 


« i 4 3,2 E) 
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2 1=0 * ! (jf'°E)) 

OO 

3 E 
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4 M) (») 
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i=0 ■'i ( M ) 


«E 


2=0 

OO 


J, 

j'P>M 

Jf 0) («)' 

(Jf’Ylf 


-T 

2 fa (J, (3 '">W) 


, 2 ’ 


The computation results of these examples are found in 
Section 


IV. Asymptotic analysis of extreme cases 


The form of the downlink coverage probability obtained 
in Theorem Q] seems complicated and it is hard to derive 
any qualitative property from it. Thus, in this section, we 
investigate the asymptotic property as 6 —> oo for two extreme 
cases; that is, m = 1 and m —» oo of the Nakagami-m fading. 
The case of m = 1 reduces to the Rayleigh fading and the 
asymptotic property is obtained in Q as follows. 

Proposition 2 (Theorem 2 in f\5f): The downlink coverage 
probability for m = 1 (Rayleigh fading) satisfies 


lim 9 1 / fs P(SIRg > 6) 

#—>•00 


OO a 

iMMk)) 


3 = 2 


du, 

( 21 ) 


where Y) ~ Gam(j, 1), j = 1,2,..., are mutually indepen¬ 
dent. 


Since Gam(m, 1/m) converges weakly to the Dirac mea¬ 
sure (5i with mass at the unity as m —>• oo, considering the 
case of m —> oo corresponds to ignoring the fading effects. In 
this case, we first show that the downlink coverage probability 
has the same asymptotic decay rate as the case of m = 1 and 
then derive the relation among the two asymptotic constants. 

Theorem 2: In the extreme case of m —» oo, the downlink 
coverage probability satisfies 


lim 6» 1 / /3 P(SIRi oo) 
#—>•00 



( 22 ) 


where Y) ~ Gam(j, 1), j = 1,2,..., are mutually indepen¬ 
dent. 


Furthermore, let c < ri and c(°°> denote the right-hand sides 
of (|2T1 > and < 1 22 b . respectively. Then, these asymptotic constants 
satisfy 

c ( 1} >r(l + i)c (oo) , (23) 

where T denotes Euler’s gamma function. 

Proof: Applying Hi = 1 and ||2G|| 2 ~ Yi for iGN, and 
also £(r) = r~ 2/j in ([TJ and Q, we have 

OO 

P(SIR£ oo) > 0) = P(SIRo° o) > 0, Bo = i) 

2=1 

00 V B 

= E P (? E (g) <^Y j >Y i ,j€ti). (24) 

»=i jew\{i} 3 

First, we consider the summand for i = 1 on the right-hand 
side of (l24l >. Applying the density function of Y\ ~ Exp(l) 
yields 



Yj > 0 1/0 v,j = 2,3,...) du, 


where the last equality follows from the substitution of v = 
u. Therefore, since e~ e /Py f 1 and l{ Yj >e-i/P t,} T 1 
a.s. as 9 —> oo, the monotone convergence theorem implies 




#—>oo 


3 =2 


and we obtain the right-hand side of d22l >. 


It remains to show that the summation over i = 2, 3, ... on 
the right-hand side of d24ll is o(0 _1 / /3 ) as 9 —> oo. For i > 2, 
applying the density function of Y) ~ Gam(i, 1) yields 

P ( d E (y;) 1 * <l,Y j >Y i ,jeK) 


ieN\{i} 3 ‘ 


j£N\{*} 13 


1 0 (* - !)■ 

r (8Zi)- 1//3 


du 


= E 


e~ u u i 1 


0 


(i~ 1) 


r dM )’ 


(25) 


where Z.- L = Therefore, since e “ < 1, 


(ez,) 


-1//3 


e “ u 


(<- 1 ) 


1 v f}—i/p 

r d U )<— E(V //J ). (26) 



















Here, noting that i ^ 1 and Y\ ~ Exp(l), 


E {Zi 


-i/ih _ 


= E 


( E £) 


1 \~i/0 


< E(Fi*) = i\. (27) 


Hence, applying ( l26l ) and ( f27l ) to ( 1 2rS b . we obtain 


°° ~V 6 

0 1//3 E P ( 0 E (yO < 1) Y o ^ Y i,j e n) 

*=2 j€N\{i} j 

OO 

< V = -==-» 0 as e ->• oo, 

^ 0VP - 1 

i=2 

which completes the proof of the first part in Theorem [2] 


For the second part of the theorem, since Y., j £ N, are 
mutually independent, the right-hand side of d21b reduces to 


c« = 


E ( ex py E lo g ( 1 + (y- 

\ 1 j=2 j 

poo 

> / E(e~ vf>Zl )dv, 

Jo 


du 


where Z\ = (1/^J^) and the inequality follows from 

log(l + a:) < x. Thus, substituting t = v 13 Z\, 



E(e~ vfiZl )dv 



Note in (|23T > that T(1 + l//3)'|'l as ft 1 and ft —> oo, 
and it takes the minimum value of the gamma function on the 
positive domain T( 7 min ) = 0.8856031944- • • at 1 + 1/ft = 
7 min = 1.4616321449 • • •; that is, the coefficient on the right- 
hand side of (l23l) is somewhat close to 1. 



Fig. 1. Comparison of downlink coverage probability between the GPP and 
PPP based models for /3 = 1.25 (top) and (3 = 2.0 (bottom). 


V. Numerical experiments 

The results of some numerical experiments are presented. 
The first experiment is the comparison between the GPP 
and PPP based models. In Figure [I] the downlink coverage 
probability over Nakagami-2 fading is plotted for both the 
GPP and PPP based models, where two cases ft = 1.25 and 
ft = 2.0 (that is, £(r) = r~ 2 5 and £(r) = r -4 ) are computed. 
For comparison, the corresponding results over the Rayleigh 
fading, which are obtained in (3), are also displayed there. As 
in the Rayleigh fading case, the downlink coverage probability 
for the GPP based model dominates that for the PPP based 
model in the Nakagami-m fading case. Furthermore, for both 
the GPP and PPP based models, the coverage probability seems 
asymptotically invariant in the value of m. This observation is, 
however, doubtful from the results of Lemma[T|and Theorem[2] 
Thus, this is further investigated in the next experiment. 

In the second experiment, we compare the downlink cov¬ 
erage probability in terms of m of the Nakagami-m fading. 
The coverage probability for m = 1, 2 and 3 is plotted in 
Figure [2] The asymptotic results obtained in Proposition [2] 
and Theorem [2] are also drawn in the same figure, where the 
asymptotic constant c^°° 1 for the case of m —> oo is estimated 
from 100 independent samples of (^°1 2 1 /Yj 13 ) 1 ^ in (l22l >. 


From the figure, indeed the downlink coverage probability 
seems increasing in m = 1, 2 and 3. However, as m —> oo, 
the asymptotic tail of the coverage probability is smaller than 
the others particularly for the value of ft close to 1. This 
observation agrees with the results of Lemma|2]and Theorem^ 
Anyway, a further investigation would be required concerning 
the impact of the value of m of the Nakagami-m fading. 
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